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We report the first double differential cross sections of two charged pions and kaons (e^e“ —>■ hhX) 
in electron-positron annihilation as a function of the fractional energies of the two hadrons for any 
charge and hadron combination. The dependence of these di-hadron cross sections on the topology 
(same, opposite-hemisphere or anywhere) is also studied with the help of the event shape variable 
thrust and its axis. The ratios of these di-hadron cross sections for different charges and hadron 
combinations directly shed light on the contributing fragmentation functions. For example, we find 
that the ratio of same-sign pion pairs over opposite-sign pion pairs drops toward higher fractional 
energies where disfavored fragmentation is expected to be suppressed. These di-hadron results are 
obtained from a 655 fb“^ data sample collected near the T(45') resonance with the Belle detector at 
the KEKB asymmetric-energy e''"e“ collider. Extending the previously published single-pion and 
single-kaon cross sections, single-proton (e'*'e“ —>■ pX) cross sections are extracted from a 159 fb“^ 
data sub-sample. 
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PACS numbers: 13.66.Bc,13.87.Fh,13.88.+e,14.20.Dh 


Quantum chromodynamics (QCD) is generally ac¬ 
cepted as the theory of the strong interaction. It de¬ 
scribes successfully many high-energy processes where 
the strong coupling is small and a perturbative treat¬ 
ment is applicable. However, the non-perturbative region 
of hadronic bound states such as the nucleon or the tran¬ 
sition from high-energetic partons into confined hadrons 
cannot be described quantitatively so far. Despite the 
recent progress in lattice QCD, both parton distribution 
and fragmentation functions (FF) remain quantities that 
need to be obtained experimentally. 

Fragmentation functions describe the density of final- 
state hadrons h from an initial parton q with a frac¬ 
tional energy z = Eh/ Eq &i a. certain energy scale Q. 
Single-hadron, unpolarized fragmentation is described by 
the function D/q{z,Q). The first z moment, summed 
over all final states, corresponds to energy conservation 
in the transition of the initial-state parton into the to¬ 
tal final state. Fragmentation functions cannot be di¬ 
rectly accessed but can be related to observable quanti¬ 
ties whenever hadrons appear in the final state. The most 
prominent connection to fragmentation functions can be 
found in single-hadron inclusive cross sections in electron- 
positron annihilation. This process provides very clean 
access as there are no hadrons in the initial state. These 
cross sections can then be related at leading order in the 
strong coupling, 05 , to fragmentation functions via 


da{e^e 


hX) 


dz 


oc 


E 


{Dlq{z,Q^)+Dlq{z,Q^)), 


( 1 ) 

where the scale Q = y/s is given by the center-of-mass 
(CMS) energy. Many such measurements have been per¬ 
formed for light hadrons at a range of CMS energies at 
the B factories [I|, 0 at LEP and SLC and other fa¬ 
cilities [Ml- The different energy scales can be related 
via DGLAP evolution 0 . 

Several global analyses of the e+e fragmentation data 
have been performed [Tsl - I^ that provide reasonable pre¬ 
cision on the sum of quark and anti-quark fragmenta¬ 
tion functions. The precise Belle and BaBar data al¬ 
lowed to reduce the uncertainties on the gluon frag¬ 
mentation function to pions, which only enters at the 
next-to-leading order in as- However, analysis of the 
single-hadron cross sections from e+e“ data alone can¬ 
not distinguish nor flavor-separate quark and anti-quark 
fragmentation and, in particular, favored vs. disfavored 
fragmentation. Favored (disfavored) fragmentation de¬ 
scribes the fragmentation of a parton into a hadron with 
(without) that parton flavor as valence content, such as 
u —>■ 7 r+ (m —>■ 7 r“). As the parton distribution functions 
are generally well known for up- and down-type flavors, 
semi-inclusive deep-inelastic scattering (SIDIS) [2l|, 
and hadron-collision [23 - [^ results provide leverage in 
truly global fits to obtain some flavor and charge sepa¬ 
rated information on pion and kaon fragmentation func¬ 


tions [ 2 ^. 

When selecting hadron pairs in e+e“ annihilation, the 
cross section can be expressed at leading order in as in 
terms of products of fragmentation functions [^ : 

d^tT(e+e“ —>■ hih2X) 
dzidz2 

E + D//^{z2)D'1/-{z,)) , (2) 

where it is assumed that both hadrons emerge from dif¬ 
ferent quarks and the scale dependence has been dropped 
for brevity. This assumption is strictly valid only at lead¬ 
ing order [s^ and for hadrons that are nearly back-to- 
back. In order to study its validity, events are analyzed 
here in three different topologies. When both hadrons 
are in the same hemisphere as defined by the thrust axis 
defined below, it is more likely that they emerge from 
the same parton so that a di-hadron fragmentation func¬ 
tion should describe the process. If both hadrons are in 
opposite hemispheres, the assumption of single-hadron 
fragmentation for each hadron is more likely. In a third 
sample, all hadron pairs irrespective of topology are con¬ 
sidered. The thrust axis h maximizes the thrust T 

The sum extends over all detected particles, and 
denotes the momentum of particle h in the CMS. 

If the assumption of single-hadron fragmentation holds 
as described in Eq. (ED, the cross sections are then sen¬ 
sitive to favored and disfavored fragmentation depend¬ 
ing on the charges and hadron types of the two detected 
hadrons. For pairs of oppositely charged pions, either 
both of the hadrons are produced by favored fragmenta¬ 
tion or both are produced by disfavored fragmentation; 
for same-sign pion pairs, one is produced from favored 
and one from disfavored fragmentation. Consequently, 
the cross section for same-sign pion pairs is smaller than 
that for opposite-sign pion pairs if disfavored fragmenta¬ 
tion functions are smaller, especially at high z as found 
in the global fits and expected in models. The reason for 
the different z dependence originates in the assumption 
that more quark-antiquark pairs need to be created to 
arrive at a disfavored hadron, which reduces its large-z 
possibility. 

When neglecting strange and charm fragmentation and 
assuming SU (2 )p isospin symmetry, the ratio of same- 
over opposite-sign pion pair cross sections reduces, for 
diagonal zi = Z 2 elements, to a simple expression of dis¬ 
favored and favored light-quark fragmentation functions. 
Strange and charm fragmentation dilute this simple re¬ 
lation for pions but in a global analysis all yields and all 
flavors can be treated appropriately. This general idea 
has been formulated in the context of the Collins frag¬ 
mentation function measurements in e“'"e“ and applied 
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there [1^ but has already been considered much 

earlier in Ref. |3Cll|. 

In the case that two hadrons are detected in the 
same hemisphere, their production is more likely to arise 
from the same parton so that di-hadron fragmentation 
functions (DiFF) should describe their yields theoreti¬ 
cally. The formalism for DiFFs was developed initially in 
Ref. [13 and including DGLAP evolution [11,113 sum¬ 
marized in Ref. [H j. Their polarized counterparts, some¬ 
times denoted as interference fragmentation functions, 
have been widely used in SIDIS experiments [H, HI] and 
Belle [13 to access together the quark transverse-spin 
distribution in the nucleon pi] . In this paper, the indi¬ 
vidual z dependence of the unpolarized baseline DiFFs 
is extracted. 

It should be noted that the leading-order mapping 
of single versus di-hadron fragmentation to the oppo¬ 
site versus same hemisphere assignments fails at next-to- 
leading order, where both t ype s of fragmentation must be 
considered simultaneously [SO]. Due to energy conserva¬ 
tion, the momentum of the same-hemisphere di-hadrons 
should not exceed the total initial momentum of the par- 
ton if originating from one parton only. 

The inclusive cross sections for charged di-hadrons in 
various topologies as a function of their fractional ener¬ 
gies Zi and Z 2 are extracted in this paper. To evaluate 
the role of favored and disfavored fragmentation, the ra¬ 
tios between these cross sections for various charge and 
hadron type combinations are calculated as well. The 
contributions for different topologies are compared to 
better understand single versus di-hadron fragmentation. 
Finally, the cross sections are compared to various Monte 
Carlo (MC) simulation tunes optimized for different col¬ 
lision systems and energies. 

Since the corrections are rather similar, a modified ver¬ 
sion of the di-hadron analysis code is used to extract 
single-hadron results as a comparison and cross check to 
the previously published single-hadron cross sections [3 ■ 
As a new result, the previously unpublished single-proton 
cross sections as a function of z are presented here and 
compared to the aforementioned MC tunes. 

This paper is organized as follows: after a short de¬ 
scription of the detector in section I, the raw di-hadron 
measurement is described in section II before detailing 
the various corrections necessary to arrive at the final 
cross sections, their ratios as well as topology depen¬ 
dence. The single-hadron analysis, including the new 
single-proton results, are presented and compared to MC 
tunes in section III. We conclude with a summary in sec¬ 
tion IV. 


I. BELLE DETECTOR AND DATA SELECTION 

This di-hadron and single-proton cross section mea¬ 
surements are based on data samples of 655 fb~^ and 
159 fb“^, respectively, collected with the Belle detector at 
the KEKB asymmetric-energy e+e” (3.5 CeV on 8 GeV) 


collider [13 operating at the T(4S') resonance (denoted 
as on-resonance) as well as 60 MeV below for comparison 
(denoted as continuum). 

The Belle detector is a large-solid-angle magnetic spec¬ 
trometer that consists of a silicon vertex detector (SVD), 
a 50-layer central drift chamber (CDC), an array of aero¬ 
gel threshold Cherenkov counters (ACC), a barrel-like ar¬ 
rangement of time-of-flight scintillation counters (TOF), 
and an electromagnetic calorimeter comprised of CsI(Tl) 
crystals (ECL) located inside a superconducting solenoid 
coil that provides a 1.5 T magnetic field. An iron flux- 
return located outside of the coil is instrumented to de¬ 
tect K'l mesons and to identify muons (KLM). The de¬ 
tector is described in detail elsewhere [2. Two inner 
detector configurations were used. A 2.0 cm beampipe 
with 1 mm thickness and a 3-layer SVD were used for 
the first sample of 97fb~^, while a 1.5 cm beampipe, a 
4-layer SVD and a small-cell inner drift chamber were 
used to record the remaining 558fb~^(159fb~^ for the 
single-hadron analysis) [il] . 

The primary light- and charm-quark simulations used 
in this analysis were generated with Pythia6.2 [13, em¬ 
bedded into the EvtCen [13 framework, followed by a 
GeantS [13 simulation of the detector response. The 
various MC samples were produced separately for light 
(uds) and charm quarks. In addition, we generated 
charged and neutral B meson pairs from T(45') decays 
in EvtCen, r pair events with the KKMC [13 generator 
and the Tauola [13 decay package, and other events 
with either Pythia or dedicated generators [^ . 


A. Event and track selection 

Events with at least three reconstructed charged tracks 
must have a visible energy of charged tracks and neutral 
clusters above 7 GeV (to remove r pair events) and either 
a heavy jet mass (the greater of the two invariant masses 
of all particles in a hemisphere) above 1.8 GeV/c^ or a 
ratio of the heavy jet mass to visible energy above 0.25. 

Tracks must be within 4 cm (2cm) of the event vertex 
along (perpendicular to) the positron beam axis. Each 
must have at least three SVD hits and fall within the 
barrel and full particle-identification (PID) polar-angle 
acceptance of —0.511 < cos^iab < 0.842. The fractional 
energy of each track must exceed 0.1. (Note that, in 
this paper, we study fragmentation functions for z above 
0.2). This initial fractional energy selection always takes 
the nominal hadron mass as given by the PID informa¬ 
tion into account. The requirement of z > 0.1 therefore 
safely accommodates pion-kaon misidentification, which 
is unfolded in the course of this analysis. 

In addition, in order to study whether two hadrons 
have likely emerged from the same parton or different 
partons, the analysis is performed on several different 
sets by requiring that both hadrons be in opposite hemi¬ 
spheres, the same hemisphere, or anywhere as depicted 
in Figs. [T] and [H For the data sets where a hemisphere 
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FIG. 1. Illustration of di-hadron fragmentation where the 
final-state hadrons are depicted as red arrows, the incoming 
leptons as blue arrows, and the event plane - spanned by 
leptons and thrust axis - is depicted as a light blue plane. 
In this case, both hadrons are found in opposite hemispheres 
defined by the thrust axis, and generally out of the plane, as 
indicated by the cones. 


assignment is required, the hemispheres are defined by 
the plane perpendicular to the thrust axis and the thrust 
must satisfy T > 0.8. 



FIG. 2. Illustration of di-hadron fragmentation where the 
final-state hadrons are depicted as red arrows, the incoming 
leptons as blue arrows, and the event plane - spanned by 
leptons and initial quarks/thrust axis - is depicted as a light 
blue plane. In this case, both hadrons are found in the same 
hemisphere as dehned by the thrust axis, and generally out 
of the plane, as indicated by the cones. 


B. PID selection 

To apply the PID correction according to the PID ef¬ 
ficiency matrices described in Ref. [l|, the same selec¬ 
tion criteria must be applied to define a charged track 


as a pion, kaon, proton, electron or muon. The informa¬ 
tion is determined from normalized likelihood ratios that 
are constructed from various detector responses. If the 
muon-hadron likelihood ratio is above 0.9, the track is 
identified as a muon. Otherwise, if the electron-hadron 
likelihood ratio is above 0.85, the track is identified as an 
electron. If neither of these applies, the track is identi¬ 
fied as a kaon by a kaon-pion likelihood ratio above 0.6 
and a kaon-proton likelihood ratio above 0.2. Pions are 
identified with the kaon-pion likelihood ratio below 0.6 
and a pion-proton ratio above 0.2. Finally, protons are 
identified with the inverse proton ratios above with kaon- 
proton and pion-proton ratios below 0.2. While neither 
muons nor electrons are considered explicitly for the sin¬ 
gle and di-hadron analysis, they are retained as necessary 
contributors for the PID correction, wherein a certain 
fraction enter the pion, kaon and proton samples under 
study. 


II. DI-HADRON ANALYSIS 

In the following sections, the di-hadron yields are ex¬ 
tracted and, successively, the various corrections and 
the corresponding systematic uncertainties are applied 
to arrive at the di-hadron differential cross sections 
(Pa{e~^e~ —>■ hih 2 X) / dzidz^- 


A. Binning and cross section extraction 

For the di-hadron cross sections, a (zi, Z 2 ) binning is 
used. We forgo a combined z and invariant-mass binning 
of the hadron pair; the latter, in particular, is relevant 
in the some-hemisphere topology as an unpolarized base¬ 
line to the previously extracted interference fragmenta¬ 
tion functions and would have allowed the extraction 
of individual fragmentation functions for p, K* ^ (j) and 
other resonances. 

The zi and Z 2 ranges of 0.2 to 1.0 used in this analysis 
are each partitioned into 16 equidistant bins. All hadron 
and charge combinations are treated independently and 
are merged only after all corrections are applied and af¬ 
ter confirming their consistency where applicable (be., 
where the same combinations of fragmentation functions 
appear, such as tt+tt''" and 7r“7r“). This leaves 16 dif¬ 
ferent charge and type combinations for pions and kaons 
initially, of which six contain irreducible information. 

Furthermore, as mentioned in the introduction, three 
hemisphere combinations are studied: two hadrons in the 
same hemisphere, two hadrons in opposite hemispheres, 
and two hadrons irrespective of hemisphere or thrust cut; 
these are abbreviated hereinafter as same, opposite, and 
any, respectively. 
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FIG. 3. (color online) Ratio of yields after and before applying the PID correction for various hadron combinations in any 
topology. For brevity, only diagonal {zi = Z 2 ) entries in each two-dimensional matrix are shown. Empty bins are visible where 
the yields become zero, especially for high-z bins. 


B. PID correction 

As in Ref. [ij, the particle misidentification is cor¬ 
rected via inverted 5x5 particle-misidentification ma¬ 
trices for the five particle hypotheses (pions, kaons, pro¬ 
tons, muons, and electrons) for each identified particle, 
laboratory momentum, and polar angle bin. These ma¬ 
trices are obtained using decays of £)*+, A and J/ip from 
data where the true particle type is determined by the 
charge reconstruction and the invariant mass distribu¬ 
tion. Occasionally, when too few events are available 
in the data, the extracted efficiencies are interpolated 
and/or extrapolated based on the behavior in the generic 
MC; this occurs particularly at the boundaries of the ac¬ 
ceptance. The matrices are calculated for each of the 
two-dimensional bins in laboratory momentum and po¬ 
lar angle, with the boundaries of the 17 bins in mo¬ 
mentum at (0.5,0.65,0.8,1.0,1.2,...., 3.0,3.5,4.0,5.0,8.0) 
GeV/c and the boundaries of the 9 bins in cos 9 
at (-0.511, -0.3, -0.152,0.017,0.209,0.355,0.435,0.541, 
0.692,0.842). 

In this analysis, the inverted misidentification matrix is 
applied for each of the identified hadrons by multiplying 
the respective weights for each hadron being a pion or 
kaon to obtain the total weight for the di-hadron and 
any of the four pion-kaon combinations. To confirm the 
consistency of this treatment, the branching ratios for 
the pion-pion and kaon-kaon decay channels to the pion- 
kaon decay channel are compared to the PDG [5^ values 


and found to be consistent. We confirm that the total 
yield of particle pairs is unaffected by this treatment. 

The corrected yields are distributed among the (zi, Z 2 ) 
bins according to the corresponding hadron masses: one 
identified hadron pair appears in several z bins with the 
above-determined weights, depending on the particular 
hadron combination. The ratios relative to the uncor¬ 
rected hadron assignment are displayed in Fig. |31 where 
one can see that the overall corrections are of the order 
of 20% to 50%. 


1. Uncertainties from the PID correction 

The uncertainties on the PID matrices are taken into 
account as uncertainties in the di-hadron yields and prop¬ 
agated through the subsequent corrections. At present, 
the uncertainties are only assigned individually for each 
hadron combination, neglecting the correlations between 
different hadron combinations. They follow the uncer¬ 
tainties assigned in Ref. [I| but take into account the 
additional complication of having two rather than one 
hadrons to unfold. To obtain the final uncertainties, 
the asymmetric uncertainties on the inverted PID matri¬ 
ces are sampled N times with a random generator with 
Gaussian distributions around the central value (sepa¬ 
rately above and below this value). From the resulting 
yields, the 68th percentiles from N samples relative to 
the central values for each (zi,Z 2 ) bin are taken as the 
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systematic uncertainties on the PID-corrected di-hadron 
yields due to the PID matrix evaluation and inversion un¬ 
certainties. The size of the statistical and systematic un¬ 
certainties relative to the PID-corrected di-hadron yields 
is displayed in Fig. |4] for selected hadron combinations. 
As expected from the overall size of the yields, the sta¬ 
tistical precision is best for opposite-sign pion pairs, fol¬ 
lowed by pion-kaon and then kaon pairs, with the same- 
sign precision being generally lower. The PID systematic 
uncertainties are mostly smaller than the statistical un¬ 
certainties in almost any bin, with the exception of the 
lowest z bins. 


C. Momentum smearing correction 

The reconstructed fractional momentum z of each 
hadron may have been smeared from its actual value and 
therefore must be corrected. For this purpose, the generic 
MC simulations are used to create two-dimensional his¬ 
tograms with 16 X 16 bins of generated and reconstructed 
izi,Z 2 ) combinations for each of the two hadrons. Only 
events that are generated and reconstructed within this 
range of z are considered. Events outside this range are 
treated in the manner described later in the acceptance 
correction section. The two-dimensional response his¬ 
tograms are created for each hadron and charge combina¬ 
tion and for all topology assignments. As the PID correc¬ 
tion was already applied to the data before the smearing 
correction, the true particle-type information in the MC 
is selected for both generated and reconstructed (zi, Z 2 ). 
In the smearing matrices, the diagonal elements are dom¬ 
inant for all hadron combinations, as can be expected in 
such a coarse binning. In the case of the pion-kaon and 
kaon-kaon combinations, the non-diagonal elements are 
slightly larger than for the pion-pion case, which indi¬ 
cates that the kaon smearing is slightly larger than that 
of pions. The small off-diagonal components facilitate the 
inversion of the smearing matrix significantly and a sim¬ 
ple, analytically inverted matrix should be sufficient to 
unfold the di-hadron yields. However, the singular value 
decomposition (SVD) unfolding method [1^ is used as 
our default to properly unfold the statistical uncertain¬ 
ties and assign systematic uncertainties due to the limited 
MC statistics, especially for bins far from the diagonal. 
This also takes into account the possible effects of the 
different z distributions in data and MC, which will be 
discussed later. 

In the SVD unfolding, a regularization parameter k 
accounts for the lack of statistics in the smearing ma¬ 
trix entries and the shape of the MC spectra. If k is 
not selected properly, the unfolded yields can be either 
too biased by the MC spectrum (too small a value for 
k) or can exhibit large fluctuations (too large a value for 
k). Following the procedure of Ref. [H^, the best reg¬ 
ularization parameter is chosen when the index of the 
regularized vector becomes smaller than unity. Initially, 
the fluctuations observed in the regularization parameter 


distribution are large, so that the false choice of too-small 
k values leads to spurious discrepancies at very high z 
between distributions of the same physics content that 
were consistent before unfolding. After smoothing the k 
distributions, the expected behavior is markedly better 
{i.e., a relevant exponentially falling contribution and an 
irrelevant flat contribution due to MC statistical fluctu¬ 
ations) and the selection of the regularization parameter 
is considered more reliable. 

The smearing is corrected in all data samples after the 
PID correction is applied and before the non-qq events are 
removed. The final before/after ratio plots are displayed 
in Fig. [S] as a function of (zi, Z 2 ). Apart from the highest 
(zi,Z 2 ) bins, where the corrections get large, the yield 
ratios are close to unity. 

All uncertainties prior to the smearing-unfolding (PID 
and statistical uncertainties) are unfolded as well, result¬ 
ing in the respective covariance matrices. The covariance 
matrix due to the MC statistics itself and the differences 
with an analytic unfolding {i.e., application of the in¬ 
verted response matrix) are assigned as systematic un¬ 
certainties related to the unfolding. 


D. Non-gg background correction 


Various QED processes can produce hadronic final 
states that contribute to our di-hadron yields and 
must be removed. Apart from particle misidentifi- 


cation, which has been addressed already, e+e" —>■ 
e+e^e+e”, and Bhabha scattering 

processes cannot contribute to hadronic final states, as 
has been verified in MC simulations. The processes that 
do produce hadron pairs are either QED processes having 
partons created initially, such as two-photon processes 
e“*"e“ — >■ e'^e~uu, e'^e~dd, e'^e~ss and e'^e~cc, or via de¬ 
cays such as from e“''e“ —^ r+T“. Hadrons from these 
processes are not produced directly via e+e" —)> qq and 
hence should not be included in our extracted di-hadron 
cross sections. Similarly, resonant T(4S') production and 
subsequent decays into neutral or charged B meson pairs 
create pion and kaon pairs that must be removed (the 
bb process does not contribute 


non-resonant e“'"e 


+ e- 


54jb The direct production of quark-antiquark pairs in 
electron-positron annihilation e“'"e“ —>■ uu, dd, ss and cc 
is treated as signal in this section, while weak decays in 
these continuum processes will be treated later. 

Figure [5] shows the relative fractions of all these pro¬ 
cesses for selected hadron pairs in the any di-hadron 
topology. Due to the large branching fraction of the 
single-prong r decay for at least one of the r leptons, r 
processes are the dominant background for pions from 
small to especially large fractional energies where the 
single hadron inherits a large fraction of the r mo¬ 
mentum. The other non-gg processes generally play a 
minor role with contributions less than a few percent, 
with the exception at high z where the two-photon pro¬ 
cess ee{uu, dd) contributes several percent to pion pairs. 




FIG. 4. (color online) Statistical (full lines), upper (dashed lines) and lower (red, dotted lines) systematic uncertainties relative 
to the PID corrected di-hadron yields in any topology as a function of {zi = Z 2 ); the diagonal bins of each two-dimensional 
matrix is shown for brevity. Empty bins are visible where the yields become zero, especially for high 2 : bins. 


Resonant T(45') production either in charged or neu¬ 
tral (CP-mixed) P meson decays contributes a few to 
about 10% and it vanishes when one fractional energy ap¬ 
proaches 0.5 due to the additional decays needed to pro¬ 
duce pions and kaons. The distributions for other di-pion 
topologies are similar except that the additional thrust 
requirement removes nearly all T(4S') decays. For same- 
sign di-pions, the r contribution is substantially smaller 
as the single-prong decays of oppositely charged r cre¬ 
ate predominantly oppositely-charged pions. For same- 
hemisphere di-pions, the single-prong r decays cannot 
contribute and consequently the relative r contributions 
are below 10-20% everywhere. 

For kaon-related di-hadron combinations, the overall 
non-qq contributions are as small as for di-pions, but 
eess and eecc are more important. In addition, the r 
decays do not play a substantial role due to the sup¬ 
pressed kaonic decays. Charm decays generally produce 
more CKM-favored kaons than CKM-suppressed pi¬ 
ons. This results in a generally larger fraction of charm 
events contributing to the pion-kaon and kaon-kaon cross 
sections: up to 60% for kaon pairs at the lowest z, with a 
similar fall-off as for pion pairs. Similarly, T(4S') decays 
favor kaons over pions and thus their fractions are as high 
as 20% (summed), rapidly disappearing at higher z. 

Assuming that the non-qq and T MC reliably describe 
the data, the background contributions are directly sub¬ 
tracted from the data distributions. In this way, we avoid 
introducing further uncertainties due to the shape of the 


udsc MC. As all these processes are QED and T(4S') pro¬ 
cesses, they are very well understood at the theory level. 
The yiel d uncertainty is 1.4% for the e+e“ —>■ t+t“ pro¬ 
cess but is substantially larger for the two-photon 
processes due to associated production, which is not 
taken into account in the current two-photon simulations. 
A factor of four relative to the nominal yield has been 
assumed for the latter [b^ . For the systematic uncer¬ 
tainties due to the non-qq background correction, these 
overall uncertainties as well as the statistical uncertain¬ 
ties in the non-qq MC, are taken into account. 


E. Preselection and acceptance correction 

The preselection and acceptance correction is divided 
into three separate terms, motivated by the different 
sources of corrections and to better expose their indi¬ 
vidual effects. The first takes into account the effect on 
the reconstruction within the specified acceptance selec¬ 
tion, mostly due to the preselection criteria and decays 
in flight; the second treats the losses outside the barrel 
acceptance; and the third takes into account potential 
losses as | cos0| approaches unity, which are not properly 
described in the generic MC. 





























9 


g 

■5 

2 

g 

■>s 


g 

■5 

2 

g 

■>s 



0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 0.3 0.4 0.5 0.6 0.7 0.8 0.9 0.3 0.4 0.5 0.6 0.7 0.8 0.9 


Z„ Zj 




Z„ Zj 


FIG. 5. (color online) Ratio of yields after and before applying the smearing correction for various hadron combinations without 
hemisphere assignment. For brevity, only diagonal {zi = Z 2 ) entries in each two-dimensional matrix are shown. 


1. Reconstruction efficiency within the barrel acceptance 


The first correction incorporates generated hadron 
pairs within the barrel geometry that do not get recon¬ 
structed. As particle identification, non-qq removal and 
smearing have already been applied, the reconstructed 
events are considered based on the generic reconstructed 
udsc MC information, but taking the MC-truth parti¬ 
cle type and momenta instead of the reconstructed val¬ 
ues. The correction factor is calculated as the ratio of re¬ 
constructed to generated events per {zi,Z 2 ) bin for each 
hemisphere assignments, hadron types, and charge com¬ 
binations. This correction takes also into account the 
events that were initially smeared out of or into the z 
range considered for this analysis. 


The efficiencies are relatively flat at around 70% and 
only drop substantially at hi^er z. This is similar to 
the behavior noticed in Ref. [J], where it was found to 
be mostly due to the preselection criteria, especially the 
heavy jet mass restriction that disfavors high-z events 
where the hadrons naturally have to be more aligned with 
the thrust axis as little other energy remains. Also, the 
minimum track requirement of three disfavors very high z 
hadron pairs, where for the same reason the multiplicity 
is small. 


2. Acceptance outside the barrel region 

A certain fraction of di-hadrons are not reconstructed 
because at least one of the hadrons is outside of the barrel 
acceptance. This fraction is evaluated by comparing the 
generated MC within the barrel acceptance {i.e., includ¬ 
ing the acceptance selection criteria) with the generated 
MC without the acceptance requirement. This accep¬ 
tance fraction is around 70% and rather flat as a func¬ 
tion of (zi, Z2), increasing slightly at very high fractional 
energies. 

The only systematic uncertainties related to these ac¬ 
ceptance corrections (both within and outside the barrel 
region) originate from the statistical uncertainty of the 
fractions within the acceptance. These uncertainties are 
rather moderate in comparison to all other systematic 
and statistical uncertainties that are aggregated in sec¬ 
tion ITTH] 


3. Large |cos0| region 

The generated MC does not necessarily reproduce 
the hadron distributions well for very forward or back¬ 
ward polar angles in the CMS. Ideally, the hadron polar 
angular distributions should resemble those of the ini¬ 
tially produced quark-antiquark pairs and thus follow a 
(1 -I- cos^ 6) dependence, neglecting the small linear de¬ 
pendence due to j — Z interference. The lower the frac- 
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FIG. 6. (color online) Fraction of hadron pairs in any topology as a function of («i, 22 ) originating from various sub-processes. 
The individual relative contributions are displayed from top to bottom for uds (red filled area), charm (blue, dotted area), 
mixed (T(4S') —>■ B^B , dark-green, hatched area) and charged (T(45') —>■ B'^B~, violet, horizontally hatched area), r pair 
(light green, scaled area), eeuu (purple, starred area), eess (light blue, dotted area) and eecc (orange hatched area) events. 
Also, for comparison, the continuum (green, solid lines) and on-resonance (orange, dotted lines) data are shown. For brevity, 
only diagonal (21 = 22 ) entries in each two-dimensional matrix are shown. 


tional energy, the less pronounced this behavior: this is 
due to the additional smearing by the transverse momen¬ 
tum generated in the fragmentation process. While such 
a behavior is roughly visible at smaller polar angles, the 
distributions rapidly drop off at higher polar angles as if 
some remaining acceptance cut is still present. As a con¬ 
sequence, the previous acceptance and efficiency correc¬ 
tions are not complete and need to be further corrected 
for this effect. As the dependence at smaller polar angles 
is well described by the expected parabola, this func¬ 
tion is used to fit the MC and compare the areas below 
the fit-result curve and the actual histograms. In prin¬ 
ciple, this treatment should be independent for the two 
hadrons and can be applied by multiplication of the two 
individual correction factors. An expected increase of the 
correction with increasing Z 2 , due to higher-z tracks be¬ 
ing more collimated and thus closer to the partonic polar 
angular dependence, has been confirmed. The overall ef¬ 
fect of this last acceptance correction is on the order of 
a few percent. 

The effect of all three acceptance and efficiency cor¬ 
rections is summarized in Fig. 0 where the ratios of the 
di-hadron yields before and after the corrections are dis¬ 
played. The overall effect amounts to between two times 


the initial yields at moderate (21,22) and more than ten 
times at very high 2 (where the event preselection cor¬ 
rection dominates). 


F. Weak decays 

Generally, fragmentation functions are only defined for 
hadrons produced by QCD processes and decays and so 
any weak decays should be removed. In practice this is 
only possible —if at all— with the help of MC and not 
entirely reliable. Therefore, many fragmentation results 
do not exclude weak decays or only those experimentally 
detectable such as those of A baryons and neutral kaons. 
The approach taken here is to provide results that either 
contain all weak decays or completely remove them with 
the help of MC. Every cc event undergoes at least one 
weak decay to produce a pion or kaon. However, in the 
fragmentation process, various quark-antiquark pairs are 
created and consequently pions and kaons can be created 
that did not originate directly from the decays of charmed 
hadrons. The only way to separate them is by following 
the parents of each final state hadron in the MC to either 
a gluonic string, which corresponds to the absence of a 
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FIG. 7. (color online) Yield ratios after successively applying all acceptance and efficiency corrections (labeled AccI to AccIII as 
the three acceptance corrections discussed in the text) relative to the reference yields before {non-qq removal) for various hadron 
combinations without hemisphere assignment. For brevity, only diagonal {zi= Z 2 ) entries in the two-dimensional matrices are 
shown. 


weak decay, or a hadron with a different, non-light va¬ 
lence flavor. In the latter case, a weak decay was present 
and this hadron would have to be removed. The difficulty 
is rapidly (algorithmically) determining this information 
for a given hadron type. In the di-hadron analysis, it 
can be argued that the chance of at least one of the 
two hadrons being from a weak decay is much higher 
for charm events and that removing all charm events is 
a valid approximation. However, this needs to be tested. 

The MC history of each hadron is studied to find weak 
decays. The heaviest flavor of each particle in the decay 
chain is selected and compared to the mother particles. 
If the decay chain ends at a string without a change in 
its heaviest flavor, no weak decay is present. If the fla¬ 
vor does change, a few strong decays need to be vetoed 
before asserting the presence of a weak decay. Examples 
are various vector mesons and other excited states with 
nonzero strangeness where the strangeness is retained in 
a lower mass state, such as K* —?> Ktt. Also, various ss 
and cc resonances need to be excluded as they also decay 
strongly despite the Zweig rule [s^ . 

The overall weak- and strong-decay fractions as a func¬ 
tion of {zi, Z 2 ) are shown in Fig. |S]for the main particle 
combinations within the any topology. Similar results are 
obtained for the other two topologies. It should be noted 
that the assumption of charm events creating only weak- 
decay pions and kaons is almost fulfilled in the procedure 


mentioned above but that a small fraction of charm di¬ 
hadron events nevertheless originates in strong decays. 
Overall, the fraction of strong decays dominates in all 
(zi,Z 2 ) bins for pion pairs, while the higher fraction of 
charm events results in larger weak fractions for pion- 
kaon and is even more pronounced in kaon-kaon combi¬ 
nations. In all cases, the weak fractions drop with 2 : as 
the additional decays soften the spectrum. 


These strong/weak fractions are mo del-dependent 
statements as the fragmentation process is only approxi¬ 
mated in Pythia. Furthermore, the absolute size of weak 
decays within Pythia depends also on the fragmentation 
settings. The uncertainties due to these effects are eval¬ 
uated by comparing the generic Belle MC to the Pythia 
default settings in the MC. The strong fractions for the 
Pythia default MC are given in the plots as well for 
comparison. As can be seen, they are rather similar but 
are generally slightly lower at high 2 : where the generally 
harder spectra in the default Pythia settings allow for 
slightly more weak decays to be present. The differences 
are assigned as a systematic uncertainty for the cross sec¬ 
tions that have the weak processes removed. 
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FIG. 8. (color online) Fraction of hadron pairs in the any topology as a function of ( 21 , 22 ) originating from weak and strong 
decays. The individual relative contributions are displayed from top to bottom for strong uds decays (purple, dark filled area), 
weak charm decays (blue dotted area), strong charm decays (dark-green, negative hatched area) and weak uds decays (red, 
horizontal striped area). The strong decay fractions are also displayed as dashed magenta lines. Also the weak decay fractions 
for udsc MC using the default Pythia settings are indicated by the dark-green, solid lines. For brevity, only diagonal (21 = 
22) entries in each of the two-dimensional matrices are shown. 


G. ISR correction 


Initial-state radiation reduces the CMS energy of the 
produced quark-antiquark pair. Consequently, the frac¬ 
tional energies calculated relative to the nominal CMS 
energy are not correct. This can alter the shape of the 
actual z dependence of the fragmentation functions and 
also invalidates pQCD calculations evaluated at the nom¬ 
inal CMS energy. The correction procedure relies on the 
strategy applied in Ref. [I| for the single-hadron cross 
sections. The events are classified according to their dif¬ 
ference from the nominal CMS energy; the events with 
a CMS energy below 99.5% of the nominal energy are 
removed. Ideally, one would want to observe the initial- 
state radiation directly in the reconstructed data; how¬ 
ever, most photons are in the very forward and backward 
regions, outside the Belle acceptance. Instead, generated 
MC data are used to directly identify ISR photons and 
remove their energies from the total CMS energy. In 
the MC, the ISR photons are identified by having their 
mother particles be an initial-state lepton. The fraction 
of such events depends on the fractional energies of the 
two final-state hadrons. If a large amount of the en¬ 
ergy is removed by the photons from the produced quark- 


antiquark system, very high fractional energies with re¬ 
spect to the nominal ^ are inaccessible. Therefore the 
fraction of non-ISR events {i.e., less than 0.5% CMS en¬ 
ergy loss) increases with increasing fractional energies. 
This is indeed the case as can be seen in Fig. 0 for the 
any topology hadron pairs (and similarly for the other 
two topologies). The events are then corrected by this 
fraction to obtain the ISR-free differential cross sections 
at the nominal center-of-mass energy. Since the ISR frac¬ 
tion depends on the fractional energies of the hadrons, 
the ( 21 , 22 ) distribution shape of the MC simulation en¬ 
ters in the ISR correction. To address the dependence of 
the ISR correction on the shape in the MC, an alternative 
MC is used for comparison and the differences in the ex¬ 
tracted di-hadron cross sections is assigned as systematic 
uncertainties. These fractions are also shown in Fig. |9l 
The ISR fractions are found to be consistent within the 
limited precision for both Pythia settings. 

The total impact of all corrections from the particle 
misidentification to the correction for ISR can be seen 
successively in Fig. [TU] for the main hadron combinations 
without a hemisphere assignment. The overall correction 
to the raw yields is substantial, predominantly due to the 
necessary acceptance corrections. They are comparable 
for most ( 21 , 22 ) bins but rise at the highest 2 bins due 
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FIG. 9. (color online) Relative Fractions of hadron pairs in the any topology as a function of ( 21 , 22 ) originating from ISR or 
non ISR events. The individual relative contributions are displayed from top to bottom for non ISR events (energy loss less 
than 0.5%, filled purple area) and ISR events (blue, dotted area) from generated generic udsc MC. The non-ISR fraction is also 
shown (red, dashed line). The relative ISR fraction for the default Pythia MC is also shown (orange, solid lines). For brevity, 
only diagonal (21 = 22 ) entries in each of the two-dimensional matrices are shown. 


to the acceptance and smearing corrections. 


H. Consistency checks and total systematic 
uncertainties 

To confirm the consistency of the results, various tests 
are performed. For example, the dependence on the data- 
taking periods is studied; after taking into account vari¬ 
ations in acceptance and reconstruction efficiency, the 
cross sections are consistent within several percent be¬ 
tween different periods and no additional systematic un¬ 
certainty is assigned. In another study, we compare the 
data recorded at the T(4S') resonance with the smaller 
off-resonance data sample. After removal of the T(4iS') 
decay contributions in the non-qq correction, the results 
from both collision energies are consistent. In yet an¬ 
other set of comparisons with the same physics-related 
information, such as charge conjugation of both parti¬ 
cles (7r+7r+ o 7T~7r~, etc.) or (random) hemisphere as¬ 
signments o K^Tr~), no systematic differences 

beyond the assigned uncertainties are found. 

All diagonal systematic uncertainties are summed in 
quadrature. The total relative systematic uncertainties 
along with the statistical uncertainties are displayed in 
Fig.[TT]for the relevant hadron pairs without topology as¬ 
signment for diagonal (zi, 22 ) bins and in Table U for the 


entire measurement range. This measurement is limited 
almost everywhere by the systematic uncertainties, for 
which the dominant contributions arise from the smear¬ 
ing correction except at high z where the rapidly falling 
MC precision contributes comparably. With increased 
MC data the systematic uncertainties could be reduced 
to the level of the statistical uncertainties and be domi¬ 
nated by the smearing correction. Additionally, there are 
global scale uncertainties due to the luminosity measure¬ 
ment (1.4%) and the track reconstruction (2x0.35%) are 
not shown. 


I. Results 

The final cross sections for the main hadron-pair com¬ 
binations are presented in Fig. [T^ for all ( 21 , 22 ) bins 
and without topology assignment. The results shown 
here and elsewhere include weak decays unless other¬ 
wise noted. As expected, the opposite-sign pion pairs 
have the largest cross sections at all 2 combinations, fol¬ 
lowed by the same-sign pion pairs. However, the op¬ 
positely charged pion-kaon and kaon-kaon combinations 
seem to be of similar magnitude or even larger than the 
same-sign pions at higher 2 , which might be explained 
by the potentially larger favored fragmentation combi¬ 
nation from strange-quark pairs. Same-sign kaon pairs 
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FIG. 10. (color online) Ratio of yields after various corrections relative to the raw yields for the main hadron combinations 
without hemisphere assignment. From top to bottom, the ratios after acceptance correction (purple, dash-dotted line), initial 
state radiation (red, long dash-dotted line), PID correction (dark green, full lines), smearing correction (blue, dashed line) and 
non-qq removal (magenta, dotted lines) are shown. For brevity, only diagonal {zi = Z 2 ) entries in each of the two-dimensional 
matrices are shown. 


TABLE I. Systematic and statistical uncertainty contributions for the main hadron combinations in the any topology inte¬ 
grated over the entire {zi,Z 2 ) range. The uncertainties due to the luminosity and track reconstruction are additional global 
uncertainties. 



1 

+ 


TT-’-yr'^ 


7V'^ 

K- 


■K+K+ 


K+K- 


K+K+ 

Statistical 

8.71 

10’ 

-05 

1.11 

10- 

-04 

1.56 

10- 

-04 

1.73 

10- 

-04 

1.83 

10- 

-04 

3.31 

10- 

-04 

PID 

9.61 

10- 

-04 

4.78 

10- 

-04 

2.09 

10- 

-03 

1.85 

10- 

-03 

2.57 

10- 

-03 

3.06 

10- 

-03 

Smearing 

6.31 

10- 

-05 

3.42 

10- 

-05 

3.92 

10- 

-04 

2.07 

10- 

-05 

6.69 

10- 

-05 

2.75 

10- 

-04 

non-qq 

6.07 

10- 

-04 

6.30 

10- 

-04 

1.03 

10- 

-03 

9.98 

10- 

-04 

1.14 

10- 

-03 

1.88 

10- 

-03 

Acceptance 

1.16 

10- 

-03 

1.32 

10- 

-03 

2.04 

10- 

-03 

2.14 

10- 

-03 

2.24 

10- 

-03 

3.65 

10- 

-03 

ISR 

3.66 

10- 

-04 

4.13 

10- 

-04 

5.97 

10- 

-04 

6.09 

10- 

-04 

7.12 

10- 

-04 

1.03 

10- 

-03 

Combined systematics 1.86 

10- 

-03 

1.71 

10- 

-03 

3.82 

10- 

-03 

4.38 

10- 

-03 

4.21 

10- 

-03 

5.28 

10- 

-02 


Luminosity 1.4 • 10 


Track reconstruction 0.7 ■ 10 


have the lowest cross sections in general, with the rela¬ 
tive differences from the other combinations increasing at 
increasing z. As at least one kaon in this case needs to be 
produced from disfavored fragmentation, the additional 
strangeness suppresses the cross sections beyond that for 
the disfavored pion fragmentation functions. 

The cross sections for di-hadrons in the same hemi¬ 
sphere are displayed in Fig. [T3] for all (^ 1 ,^ 2 ) bins. The 


cross sections fall off rapidly and mostly disappear at the 
boundary Z 1 +Z 2 = 1, where the total energy of one initial 
parton is fully contained in the energy of the two hadrons. 
A small excess above this limit can be seen. MC studies 
show that this excess can be explained qualitatively by a 
small mis-assignment of hemisphere due to the smearing 
of the thrust axis relative to the initial quark-antiquark 
axis. In addition, hard gluon radiation may create such 
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FIG. 11. (color online) Relative (asymmetric) systematic (lower uncertainties - dashed lines; upper uncertainties - dash-dotted 
lines) and statistical nncertainties (full lines) for the most relevant hadron pairs in the any topology as a fnnction of (zi,Z2)- 
For brevity, only the diagonal bins (zi = Z2) are shown. 


events. 


1. Cross section ratios 

As several of the uncertainties are common to all 
charge and hadron combinations, these cancel in ratios 
and the information about favored and disfavored frag¬ 
mentation should be more reliable. For example, ignoring 
strange-quark fragmentation to pions, the same-sign to 
opposite-sign pion-pair ratios are simple measures of dis¬ 
favored vs. favored pion fragmentation functions for light 
quarks. As can be seen in the ratios in Fig. 1141 they show 
a similar nearly flat behavior at low fractional energies, 
but deviate substantially from this trend at higher frac¬ 
tional energies. Generally, the opposite-sign pion-kaon 
and kaon-kaon ratios are more suppressed than the al¬ 
ready disfavored same-sign pion pairs, but at the high¬ 
est (zi,Z 2 ) do become comparable. This might be re¬ 
lated to the fact that opposite-sign kaons can be favored 
fragmentation for both hadrons if a ss was created ini¬ 
tially. As the opposite-sign pion-kaon pairs are similar 
in magnitude, even at high 0 , it can be argued that dis¬ 
favored fragmentation from strange quarks to pions is 
not as suppressed. Same-sign kaon pairs, where at least 
one strange-quark pair needs to be produced in the frag¬ 
mentation, are always suppressed at least one order of 


magnitude relative to the opposite-sign pion pairs. This 
shows that strangeness produced in fragmentation is in¬ 
deed strongly suppressed, as is generally assumed in frag¬ 
mentation models such as those included in Pythia. 


2. Hemisphere deeomposition 

Figure [T5] displays all six relevant hadron combinations 
for opposite hemispheres while Fig. [16] shows the cross 
sections for hadrons within the same hemisphere using a 
thrust requirement of T > 0.8. Note that the require¬ 
ment of a minimum thrust value is not corrected for in 
these hemisphere decompositions, which must be taken 
into account when used for global FF analyses. As ex¬ 
pected, the cross sections at small 2 are of similar magni¬ 
tude between the same and opposite hemispheres, while 
at higher z only opposite-hemisphere pairs remain. 

These cross sections with hemisphere assignment can 
be compared to the cross sections without the hemisphere 
assignment and without the thrust requirement, as shown 
in Fig. |T7| for diagonal (^ 1 , 22 ) bins. As expected, the 
same hemisphere contributions are comparable to those 
from opposite hemispheres at very low z while they vanish 
rapidly with increasing z. This is in agreement with the 
assumption of the same hemisphere di-hadrons emerging 
predominantly from the common initial parton. If so, 
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FIG. 12. (color online) Differential cross sections for 7r''’7r“ (black circles), (blue squares), K~ (green triangles), 

(pnrple diamonds), K~ (red crosses) and (violet downward triangles) pairs in the any topology as a fnnction of 22 

for the indicated 21 bins. The error boxes represent the systematic uncertainties. 


the sum of the two hadron’s fractional energies cannot 
exceed unity and they do indeed drop to zero at ( 21 , 22 ) 
each around 0.5. 

The small difference between the sum of the same and 
opposite hemispheres and the any topology assignment 
at low 2 is due to the additional thrust selection used 
to identify hemispheres. The small deviations seen oc¬ 
casionally at high 2 are related to variations from the 
smearing unfolding and are consistent within the uncer¬ 
tainties that are not shown in this figure. 


3. MC generator comparison 

The di-hadron cross sections in the any topology are 
compared to various fragmentation settings within the 
PYTHiA/JetSet MC generator. These are displayed in 
Fig. [THl for the main hadron combinations and diago¬ 
nal ( 21 , 22 ) bins. The settings correspond to the de¬ 
fault Pythia, those currently used in Belle for the fully 
tracked GEANT simulations, as well as various other set¬ 
tings tuned to specific experiments, collision systems, and 
energies from the LEP/Tevatron, ALEPH, and HERMES 
environments. Generally, all parameterizations agree at 
very low 2 as the total production yields for certain parti¬ 
cles are best known. However, at very high 2 , the distri¬ 


butions differ greatly. It appears that for all six hadron 
and charge combinations the default Pythia and the lat¬ 
est Belle settings describe the data best even at large 
fractional energies. The LEP-based tunes generally over¬ 
shoot the data at high 2 while the older Belle and HER¬ 
MES tunes fall off much too rapidly. 


In addition to the any topology combination, which is 
dominated by the opposite topology at higher fractional 
energies, the same hemisphere combination is compared 
to these MC tunes. In principle, the different hemisphere 
combinations are sensitive to different parameters in the 
Pythia settings. An example of the comparison for 
Tv'^K~ pairs is displayed in Fig. 1191 other hadron com¬ 
binations are available in the supplementary file. The 
overall behavior is similar to Fig.[THl with the older Belle 
and HERMES tunes undershooting the data. The other 
parameterizations do not differ as substantially as in the 
any hemisphere combinations and they all reproduce the 
data reasonably well. The differences in the Pythia set¬ 
tings are summarized in Table 1 of the supplementary 
file. 
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FIG. 13. (color online) Differential cross sections for 7r''’7r“ (black circles), (blue squares), K~ (green triangles), 

(pnrple diamonds), K~ (red crosses) and (violet downward triangles) pairs in the same topology (including T > 0.8 

selection) as a function of 22 for the indicated zi bins. The error boxes represent the systematic uncertainties. The vertical, 
dashed green line represents the 21 + 22 = 1 limit in each panel. 


III. SINGLE-HADRON ANALYSIS 

In addition to the di-hadron analysis, the production 
of single hadrons, especially previously unpublished sin¬ 
gle protons is considered here. The single-hadron anal¬ 
ysis follows the same procedure as the di-hadron anal¬ 
ysis. The 2 range between 0.1 and 1.0 is divided into 
36 bins; for protons, 2 < 0.2 is kinematically inacces¬ 
sible. The particle misidentification correction is per¬ 
formed as in the di-hadron analysis (but only for one 
track) and the resulting yield modification is shown in 
Fig. [2ni At intermediate 2 , in particular, the proton 
yields are reduced substantially due to proton misidenti¬ 
fication. Non-gg events contribute once again to the pion 
and kaon distributions but not as much to protons, where 
predominantly eeuu processes at high 2 (~ 5%) and T 
decays at low 2 (maximally ~ 20%) are the dominant 
backgrounds. All acceptance corrections are only weakly 
dependent on hadron type and show the same moderate 
(substantial) correction factors at small and intermedi¬ 
ate (high) 2 ; the high -2 correction is again dominated by 
the event preselection efficiencies. Weak decays originate 
predominantly from charm decays and so are a very small 
contribution (< 10%) for protons. The various correction 
steps for single-pions, kaons and protons are summarized 


in Fig. [201 

The ISR correction here is similar to that in the di¬ 
hadron analysis. To clarify the correction for the previ¬ 
ous single-pion and kaon results [I|, we show in Fig. |2I] 
the ISR and non-ISR fractions for single pions and kaons 
as well as protons. As in the di-hadron analysis, the frac¬ 
tion of events with an actual CMS energy below 99.5% 
of the nominal energy is below 30% and decreases with 
increasing 2 . 

The resulting single-pion, -kaon, and -proton cross sec¬ 
tions are displayed in Fig. |22l While the pion and kaon 
results are consistent within uncertainties to those pub¬ 
lished before, the proton results from Belle are shown 
for the first time. The results are compared with the 
aforementioned PYTHiA/JetSet fragmentation tunes in 
Fig. 1231 As has been noted above and in [l|, the 
PYTHiA/JetSet settings close to the default settings re¬ 
produce the pion and kaon cross sections rather well. For 
the proton cross sections, no setting describes the en¬ 
tire 2 range, while the ALEPH and LEP/Tevatron tunes 
roughly agree with the data at low 2 and the older Belle 
MC setting is in moderate agreement at high 2 . 
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FIG. 14. (color online) Ratios of differential cross sections for (blue squares), 7r''’R'“ (green triangles), -k^K'^ (purple 

diamonds), K'^K~ (red crosses) and (violet downward triangles) pairs relative to the 7r''"7r~ pairs in the any topology 

as a function of 22 for the indicated 21 bins. The shaded areas correspond to the systematic uncertainties under the assumption 
that they are independent between each two di-hadron combinations. 


IV. SUMMARY 

In summary, we present e“'"e“ —>■ hih 2 X differential 
cross sections in 21 and 22 for pion-pion, pion-kaon and 
kaon-kaon pairs of the same and opposite charges and 
in various topologies. The general expectations of dis¬ 
favored fragmentation functions being suppressed, espe¬ 
cially at large fractional energies, are confirmed within 
the assumptions mentioned in this article. In particu¬ 
lar, the same-sign pion pairs in opposite hemispheres fall 
off more rapidly than the opposite-sign pion pairs. The 
ordering with additional strangeness is also as expected 
when taking into account the favored-kaon fragmentation 
of strange quarks and charm decays. For example, where 
strangeness needs to be created in the fragmentation such 
as for same-sign kaon pairs and, to a lesser extent, the 
same-sign pion-kaon pairs, the cross sections decrease 
even more rapidly as the already disfavored same-sign 
pion pairs. 

The vanishing of the same-hemisphere di-hadron cross 
sections once the sum of the fractional energies of the 
two hadrons exceeds unity supports the assumption of 
the same hemisphere di-hadrons being produced predom¬ 
inantly via single-parton di-hadron fragmentation. This, 
in turn, bolsters the interpretation of the opposite hemi¬ 


sphere hadron pairs as arising from the fragmentation 
of different partons. As a consequence, the inclusion of 
the opposite hemisphere di-hadrons in terms of single¬ 
hadron fragmentation into a global pQCD fragmentation 
function analysis should be possible. 

The extracted di-hadron cross sections are compared 
to various Pythia MC tunes that were optimized for var¬ 
ious other energies and collision systems. A full optimiza¬ 
tion at Belle energies should be possible based on these 
results; nevertheless, both the default Pythia fragmen¬ 
tation setting as well as the latest Belle fragmentation 
setting already describe the data reasonably well. 

Single-proton e+e“ —i pX cross sections differential in 
2 are presented in addition to the previously published 
single-pion and -kaon results. These are expected to be of 
use in global analyses of fragmentation functions, includ¬ 
ing the proton results from BaBar Q. Various Pythia 
tunes are compared with conclusions similar to those in 
the di-hadron case. However, for the proton cross sec¬ 
tions, the agreement is fair at best over the entire 2 range, 
which suggests room for improvement in the Pythia set¬ 
tings to better model the baryon production. 

With the precision of these measurements and the ad¬ 
ditional information obtained by the use of di-hadrons, 
we expect that subsequent global fits to the world data 
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FIG. 15. Differential cross sections for the main hadron pairs within opposite hemispheres (including thrust selection T > 0.8) 
as a function of {zi,Z 2 ) for diagonal bins only. The error boxes represent the systematic uncertainties. 


will improve substantially our understanding of fragmen¬ 
tation functions, in particular in terms of the distinction 
of favored versus disfavored fragmentation. 
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FIG. 16. Differential cross sections for the main hadron pairs within the same hemisphere (including thrust selection T > 0.8) 
as a function of {zi,Z 2 ) for diagonal bins only. The error boxes represent the systematic uncertainties. 
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FIG. 17. (color online) Differential cross sections (top) and ratios to the same and opposite topology sum (bottom) for the main 
hadron combinations, stacking same (gray filled areas, including thrust selection T i 0.8) and opposite (blue hatched areas, 
including thrust selection T i 0.8) hemisphere data and comparing to those without hemisphere assignment (red curves). For 
visibility, only the diagonal ( 21 , 2 : 2 ) bins and statistical uncertainties are displayed. 
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FIG. 18. (color online) Differential cross sections (top) and ratios to the data (bottom) for the main hadron pairs in the any 
topology as a function of {zi, Z2) for diagonal bins only. Various Pythia tunes are also displayed as described in the text. For 
comparison, the relative statistical and systematic uncertainties are shown for the data as well. 
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FIG. 19. (color online) Differential cross sections for K pairs in the same hemisphere as a function of ( 21 , 22 ). Various 
Pythia tunes are also displayed as described in the text. 



FIG. 20. (color online) Ratio of yields after various corrections relative to the raw yields for the main single hadrons are shown 
as a function of 2 . From top to bottom, the ratios after acceptance correction (purple, dash-dotted line), initial state radiation 
(red, long dash-dotted line), PID correction (dark green, full lines), smearing correction (blue, dashed line) and non-qq removal 
(magenta, dotted lines) are shown. 
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FIG. 21. (color online) Relative fractions for pions, kaons and protons as a function of Zi originating from ISR or non ISR 
events. The individual relative contributions are displayed from top to bottom for non ISR events (energy loss less than 0.5%, 
purple, filled area) and ISR events (blue, hatched area) from generated generic udsc MC. The relative fractions are also shown 
for the non-ISR fraction (green, dashed lines). 



FIG. 22. (color online) Single-pion (black circles), -kaon (blue 
squares) and -proton (green triangles) cross sections from top 
to bottom, as a function of 2 . 
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FIG. 23. (color online) Differential cross sections (top) and ratios to the data (bottom) for the main single-hadrons as a function 
of 2 . Various MC tunes are also displayed as described in the text. For comparison, the relative statistical and systematic 
uncertainties are shown for the data as well. 














